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Abstract. In this talk we summarize the status of theoretical predictions for the 
average number of charm quarks in a B-hadron decay. 



1. Introduction 

Since quite a long time there exists a discrepancy between theoretical predictions and 
measurements of the quantity ric, which describes the average number of charm quarks 
in the final state of a B-hadron decay [Q]. In the last years this difference became smaller 
and it became a matter of taste whether one speaks of a missing charm puzzle or not. In 
this talk we try to summarize the theoretical results and to clarify the origin of different 
numbers for Uc- 

One can calculate ric in the following ways: 

nc = + ^- ^ + 2^- (1) 

i tot J- tot 

^ ^ ^ T{h^ 2c) _ V{h^ Oc) 

r(6 ^ Oc) sums up all charmless decay rates like the non-leptonic channels h uus, d, 
b sss, d,b ^ dds, d and the semi-leptonic channels b ulv and b sg, gg. r{b Ic) 
sums up all decay rates with one charm quark in the final state, like the non-leptonic 
channels b — > cus, d, b ucs, d and the semi-leptonic channels b —>■ civ. Finally we 
have r(6 ^ 2c) with two charm quarks in the final state: b — > ccs, d. 
Before we compare experimental results and theoretical predictions, let us look at the 
calculation of these decay rates. 



2. Calculation of inclusive decay rates 



The Heavy Quark Expansion (HQE) (for a recent review see 0) is the theoretical 
framework to handle inclusive S-decays. It allows us to expand the decay rate in the 
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following way 

r = ro+f— )'r2+f— )'r3 + --- (4) 

Here we have an systematic expansion in the small parameter A/m?,. The different terms 
have the following physical interpretations: 

• Fq: The leading term is described by the decay of a free quark (parton model), we 
have no non-perturbative corrections. 

• Fi: In the derivation of eq. (|^) we make an operator product expansion. From 
dimensional reasons we do not get an operator which would contribute to this 
order in the HQE. | 

• F2: First non-perturbative corrections arise at the second order in the expansion 
due to the kinetic and the chromomagnetic operator. They can be regarded as the 
first terms in a non-relativistic expansion. 

• F3: In the third order we get the so-called weak annihilation and pauli interference 
diagrams. Here the spectator quark is included for the first time. These diagrams 
give rise to different lifetimes for different B hadrons. 

• The dots represent higher order terms in l/mt, possible non-perturbative 
corrections (like in the decay B — > Xg'-f ^) and unknown terms which are due to 
duality violation (see @] for a nice review). 

Schematically one can write the Fj's as products of perturbatively calculable functions 
(depending on couplings, masses, renormalization scale,...) and matrix elements, which 
have to be determined by some non-perturbative methods like lattice-QCD or sum rules. 
Now we may have a closer look at eq. @. Each of the appearing terms can be expanded 
in a power series in the strong coupling constant 

F, = Ff ) + ^F« + . . . . (5) 

We start with a discussion of the perturbative part of the Fp^'s and then we make some 
comments about the status of the non-perturbative parameters. 

2.1. Leading term: Fq 

Fq"^ is well known. In addition we have analytic expressions of Fo^"* for b — > clu 
and b cud ^ and a numerical value for b ccs 0]. The effects of the charm 
quark mass were found to be quite sizeable. Although suppressed by one power of 
a^, penguin diagrams are dominant for b ^ no charm P], |^. Recently the NLO 
calculation for b sg has been finished [|^. The inclusion of penguin diagrams with 



current-current operators for the decay b ccs and penguin diagrams with penguin 
operators for b no charm is still missing, but their effects are not expected to be 

J Strictly spoken we get one operator of the appropriate dimension, but with the equations of motion 
we can incorporate it in the leading term. 
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large. It is a remarkable feature of the HQE that in the leading term Fq only the unit 
operator appears, so the matrix elements of this operator are trivial. Therefore we have 
no non-perturbative parameters in Fq. 



2.2. Sub-leading term: F2 



F2'''* is known for the most important operator insertions |TT]]. Some penguin operator 



insertions are still missing. It would be nice to have a result for Fg^^ but the calculation 
seems to be quite tough. One has to calculate the imaginary part of three loop diagrams 
with one external gluon. Here we have two matrix elements: Ai and A2. The first one 
is not very well known, see e.g |12|, while the second number can be extracted from 
experiment. 



2.3. Spectator effects: F3 

Spectator effects arise first in the third order of the expansion in 1/mi,. vf'' is known for 
AFBg[0 and for Bs and with charm quark mass effects [Q. F3^'' was calculated 



for AF^s by [|T5|. The calculation of F^'^ for 5+, B, and A,, is still missing. In Fq we 
have the following non-perturbative parameters: decay constants /m (depending on the 
decaying meson M) and Bag-Barameters Boj^i (depeding on the decaying meson M and 
the Dirac structures D of the appearing operators). For AF^^ we have already quite 
stable lattice predictions for these quantities, while for 5"*", Bg and A;, relieable numbers 



are still missing (see |]16|, |T7| ). 



2.4. ^/fnt corrections: F4 

For ATbs even vf^ has been calculated by This could be done for 5+, Bg and Af,, 
too. The appearing matrix elements were estimated in vacuum insertion approximation. 



3. Different normalization 



In order to determine Uc we have to determine the branching ratios for b decays into 0,1 
and 2 charm quarks. So one could simply calculate F(0, 1, 2c) and Tfot- But there are 
several reasons, why it might be better not to calculate these quantities straightforward. 
First, the semi-leptonic decay rate Tgi is clearly the most reliable prediction, while Ttot 
is probably the least reliable prediction. By writing 

B,^x = ^*^=:rx* B^^ (6) 

i- si J- tot 

we can eliminate Ttot in favor of F^^. In rx we have no m\- and Ai-dependence anymore. 
Second, the decay b — > ccx is most sensitive to possible quark hadron duality violations. 
This is due to the fact that the HQE is actually not an expansion in 1/mb, but in 
1/ where E is the energy release in the decay. For b — ccx we have E = — 2mc, 
which is already quite a small number. If we use eq. (3) and the r's instead of the 
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branching ratios, we have ehminated the decay h — ccx, as proposed in |T^. Now r(Oc) 
is an important input parameter for the determination of Uc- Possible enhancements 
of r(Oc) due to new physics would lower B*^^°^^ and •n}^^°^y simultaneously. Different 



mechanisms for such an enhancement were studied in the literature ||2^ 



4. Results in the literature 

Now we summarize the results for the relevant decay rates from the literature and 
determine Uc in various ways. 



4.1. Counting of one Charm Quark 

The dominant decay is 6 ^ cud. There was quite a confusion due to two different 
numbers in the literature: Ball et al. quote r{cud) = 4.0 ± 0.40, while Neubert 
was showing r{cud) = 4.2 ± 0.4[^ in Jerusalem. The difference of these numbers 
is an effect of second order in as- While the authors of were calculating ratios like 
(a + asb)/{c + Usd) nummerically, the author of expanded the ratio in lH- 
Unfortunateley the difference is quite sizeable. For all possible semi-leptonic decays 
we get r^iu = 2.22 ± 0.04 and for the Cabibbo suppressed decay modes the result is 
rues' = 0.03 ± 0.00. Depending on our input for r{cud) we get two different results: 



r(lc) = 6.25 ±0.4 



r(lc) = 6.45 ±0.4 ^ 



4.2. Counting of no Charm Quark 

For the non-leptonic charmless 6-decays it turned out, that penguin diagrams are as 
important as the leading contribution to these decays, although being suppressed by 
Us [@]. Even contributions, so-called double penguins have a sizeable value 
One gets r(Oc) = 0.18 ± 0.08 g, i for all charmless final states. Recently the NLO 
QCD calculation of b sg and b —>■ sgg was finished ||lOi- Greub and Liniger get an 



enhancement of more than 100% compared to the LO value 

, f 0.022 ±0.008 LO 
<^-^-^9.sgg) = ^ 0.05 ±0.01 NLO. 

With the new result for b —>■ sg and b —>■ sgg at hand we get: 
r(Oc) = 0.21 ±0.08 PI 



4-3. Counting of two Charm Quarks 

For b —>■ CCS we have again two different results. Ball et. al quote r(2c) = 2.0 =F 0.5 



while Neubert gets 1.89 ^ 0.54 21 . The difference has the same origin as in section 4.1. 
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4-.4- Results for Uc 

With the experimental value for the semi-leptonic branching ratio presented in Osaka 
^si^ = 0.1059 ± 0.0016 we can determine Uc in three different ways. 

(i) Elimination of no charm: = (r(lc) + 2r(2c)) B^J^ = 1.09 ±0.11 

(ii) Elimination of one charm: Uc = 1 + (r(2c) — r(Oc)) B'^f^ = 1.18 ± 0.06 

(iii) Elimination of two charm: = 2 — (r(lc) + 2r(0c)) -B^f^ = 1.28 ± 0.05 

For r(lc) and r(2c) we used the average of ^ and |2l|. Of course, all these numbers 
should be the same. The reason for the disagreement is found by comparing the 
theoretical and experimental value of the semi-leptonic branching ratio. Theory tells us 

(r(Oc) + r(lc) + r(2c))-' = 0.118 ± 0.009 = . 

The central value is quite above the experimental number for Bsi, but the errors are 
large. When we introduced rx in eq. (|]), we asummed that Bl'l'^"^^ = B^f^', which is 
not satisfied. This is the reason for the inconsistencies in the determination of Uc- If we 
use B^^'^°^^ to determine nc, we get in all three cases the central value = 1.21. 
In Osaka ric = 1.16 ± 0.05 was given as the experimental value [Q, while Kagan gets a 
value of ric = 1.085 ± 0.05 [Q. It is beyond the scope of this talk to clarify the origin 
of these two different experimental numbers . 



5. Disscussion and outlook 



In this talk we tried to clarify the orgin of different values for ric on the market. First 
we have different numbers for r(lc) and r(2c) due to a different treatment of 
contributions. The numbers of |^T| give a slightly smaller value for Uc, than the numbers 
of 1^, 01 . Second, we get quite different results for the three possibilities (eq. (1)- 
(3)) to determine if we use a normalization of the decay rates to F^; instead of 
Ttot^- The reason for that is the disagreement of the theoretical number for Bgi with 
the experimental value. This problem has to be resolved in the future. Third, the 
experimental value of Uc seems to be not completely clear. 

So we are still not in the position to say the final word about the existence of a 
missing charm puzzle. If we use an appropriate theoretical input and set n = mi,/ 4, 
(which means a high value for and mc/mh = 0.33|]], than experiment (the numbers 
shown in Osaka) and theory agree more or less. On the other hand there is still room for 
a deviation, which might be due to a new physics enhanced r(Oc) or dualtity violation 
in b —>■ CCS or.... Precise experimental values of r(2c) and r(Oc) would help a lot, to 
confirm or to rule out these interesting possibilities. 

§ In the determination of Ar^^, we have the same situation, that we get quite different numbers for 
different normaUzations (see talk p^). 

II Here one should keep in mind, that the ratio rric/mb is fixed by HQET. 
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